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Abstract. This paper is concerned with the study of the wave equation on compact 
surfaces and locally distributed damping, described by 

u u — &mu + a(x) g(u t ) =0 on M x ]0, oo[ , 

where M C R 3 is a smooth (of class C 3 ) oriented embedded compact surface without 
boundary, such that M = Mo U Mi, where 

Mi := {x e M; m(x) ■ v(x) > 0} and M = M\Mi. 

Here, m(x) :— x — x°, (x° 6 R 3 fixed) and v is the exterior unit normal vector field of 
M. 

For i — 1, . . . , k, assume that there exist open subsets Moi C Mo of M with smooth 
boundary dMoi such that Moi are umbilical, or more generally, that the principal curva- 
tures fei and &2 satisfy \ki(x) — kz(x)\ < £; (e< considered small enough) for all x £ Moi- 
Moreover suppose that the mean curvature H of each Moi is non-positive (i.e. H < on 
Moi for every i — 1, . . . ,k). If a(x) > ao > on an open subset M* C M that contains 
A4\ U* =1 Moi and if g is a monotonic increasing function such that k\s\ < \g(s)\ < K\s\ 
for all \s\ > 1, then uniform decay rates of the energy hold. 

1. Introduction 

Let M be a smooth (of class C 3 ) oriented embedded compact surface without boundary 
in R 3 with M = Mo UMi, where 

(1.1) M\ := {x G M;m(x) ■ u{x) > 0} and Mo = M\M\. 

Here, m(x) := x — x°, (x° G M 3 fixed) and v is the exterior unit normal vector field of M. 

We denote by Vt the tangential-gradient on M, by A_a^ the Laplace-Beltrami operator 
on M. This paper is devoted to the study of the uniform stabilization of solutions of the 
following damped problem 

uu — A_A/jii + a(x) g{ut) = on M x ]0, oof , 

u(x, 0) = u°(x), ut(x , 0) = u 1 (x) x £ M, 

where a(x) > ao > on an open proper subset M* of M and in addition g is a monotonic 
increasing function such that k\s\ < \g(s)\ < for all |s| > 1. 
Stability for the wave equation 

(1.3) u tt -Au + f(u) + a(x) g{u t ) =0infixR + , 



(1.2) 
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where f2 is a bounded domain in M n , has been studied for long time by many authors. 
When the feedback term depends on the velocity in a linear way Zuazua [ZUA] proved 
that the energy related to the above equation decays exponentially if the damping region 
contains a neighborhood of the boundary d£l of O or, at least, contains a neighborhood 
uj of the particular part given by {x € <9$7 : (x — xq) ■ v(x) > 0}. In the same direction, 
but when / = 0, it is important to mention the work due to Rauch and Taylor [Ra-Ta| 
and, subsequently, the results of Bardos, Lebeau and Rauch [BAR] . based on microlocal 
analysis, that ensures a necessary and sufficient condition to obtain exponential decay, 
namely, the damping region satisfies the well known geometric control condition. The 
classical example of an open subset u verifying this condition is when u is a neighborhood 
of the boundary. Later, again considering / = 0, Nakao [Nal| INa2] extended the results 
of Zuazua [ZUA] treating first the case of a linear degenerate equation, and then the 
case of a nonlinear dissipation p(x,u t ) (here, again, / = was considered) assuming, as 
usually, that the function p has a polynomial growth near the origin. Martinez [Mar] 
improved the previous results mentioned above in what concerns the linear wave equation 
subject to a nonlinear dissipation p(x, Ut), avoiding the polynomial growth of the function 
p(x, s) in zero. His proof is based on the piecewise multiplier technique developed by Liu 
[Liu] combined with nonlinear integral inequalities to show that the energy of the system 
decays to zero with a precise decay rate estimate if the damping region satisfies some 
geometrical conditions. More recently, and still considering / = 0, Alabau-Boussouira 
[ALA] extended the results due to Martinez [Mar] by showing optimal decay rates of 
energy. In addition, we would like to mention the most recent work in this direction due 
to D. Toundykov [Touj which presents optimal decay rates for solutions to a semilinear 
wave equation with localized interior damping and a source term, subject to Neumann- 
type boundary condition. 

A natural question arises in the context of the wave equation on compact surfaces: It 
would be possible to stabilize the system by considering a localized feedback acting on a 
part of the surface ? In affirmative case, what would be the geometrical impositions we 
have to assume on the surface? When the damping term acts on the whole surface, the 
conjecture was studied by Cavalcanti and Domingos Cavalcanti in [CA-DQ] and also by 
Andrade et al. in [Anl] IAn2] in the context of viscoelastic problems. However, as far as 
we are concerned, there is no result in the literature regarding the nonlinear wave equation 
on compact surfaces when the damping term acts in a portion A1* strictly contained in 
A4. For the linear case, we can mention the works due to Rauch [Ra-Ta] . Hitrik [HIT] and, 
more recently Christianson [CHR]. 

The main goal of this paper is exactly to prove the above conjecture when the portion 
of A4, where the damping is effective is strategically chosen. For i = l,...,k, assume 
that there exist open subsets Aioi C M.$ of M. with smooth boundary dMoi such that 
A4oi are umbilical, or more generally, that the principal curvatures k\ and ki satisfy 
\k\[x) — k%{x)\ < Ej (ej considered small enough) for all x € Aioi- Moreover suppose 
that the mean curvature H of each Moi is non-positive (i.e. H < on A^oi for every 
i = 1, . . . , k) and that the damping is effective on an open subset A4* C A4 that contains 
M\ul =1 M oi . 
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The strategy used to prove the above conjecture is basically to make use of multipliers 
and fields as in Lions [LIONS 1] with new ingredients that will be clarified in section 
4. Indeed, the main difficulty and the novelty in these kind of problems on surfaces is 
how to deal with (or to interpret) the new terms which appear in the computations that 
come from the geometrical structure of M. Moreover, this approach can be naturally 
extended for semilinear waves where the semilinear function f(s) is assumed to be super- 
linear. We would like to emphasize that the proofs of [Ra-Taj IBAR} IHIT| . based on 
microlocal analysis, do not extend to the nonlinear problem fjl .2[> . In addition, making use 
of arguments due to Lasiecka and Tataru |LA-TA| we obtain optimal decay rates of the 
energy. The obtained decay rates are optimal, since when we are able to explicit them 
(as in Cavalcanti, Domingos Cavalcanti and Lasiecka [CA-DO-LA] ). they are the same as 
these optimal rates derived in the recent works of Alabau-Boussouira [ALAj or Toudykov 

Our paper is organized as follows. Section 2 is concerned with the statement of the 
problem and we introduce some notation . Our main result is stated in Section 3. Section 
4 is devoted to the proof of the main result. 

2. Statement of Problem 

Let M be a smooth oriented embedded compact surface without boundary in R 3 with 
M = Mq U Mi, where 

(2.1) M\ := {x G M;m(x) ■ u(x) > 0} and M = M\Mx- 

Here, m is the vector field defined by m(x) := x — x°, (x° £ M 3 fixed) and v is the exterior 
unit normal vector field of M. 

In this paper, we investigate the stability properties of functions [u(x, t), ut(x, t)} which 
solve the following damped problem: 

J u tt - A M u + a(x)g{u t ) = on Aix]0,oo[, 

( ' j \u(0)=u°, u t (0)=u\ 

where the feedback function g satisfies the following assumptions: 

Assumption 2.1. 

(i) g (s) is continuous and monotone increasing, 

(ii) g (s) s > for s / 0, 

(Hi) k \s\ < g (s) < K \s\ for |s| > 1, 

where k and K are two positive constants. 

In addition, to obtain the stabilization of problem (|2.2I) . we shall need the following 
geometrical assumption: 

Assumption 2.2. Remember that for i = 1, . . . , k, Moi C Mo are open sets with smooth 
boundary dMoi such that H < and Moi are umbilical submanifolds, or more generally, 
that the principal curvatures ki and k2 satisfy \ki(x) — k2(x)\ < £j fa considered small 
enough) for all x S Moi- We assume that a € L°°(A4) is a nonnegative function such that 

(2.3) a(x) > ao > 0, a. e. on M*, 

where M* is an open set of M that contains M\ uf =1 Moi. 
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In order to fix ideas, Figure 1 shows a compact surface M such that there exists only 
one subset Moi, which we take as the interior of .Mo- 




Figure 1 

Figure 1. The observer is at xq. The subset A4q is the "visible" part of 
M and M\ is its complement. The subset M* D M — U^ =1 Moi = M\Mq± 
is an open set that contains A4\ U* =1 M-oi and the damping is effective 
there. 



In the sequel we define by £ = M x ]0, T[ , X* = M t x ]0, T[ , i = 0, 1. 

Let us considerer the Sobolev spaces H s (A4),s £ R as in Lions and Magenes [LiMa] 
section 7.3. 

On the other hand, by using the Laplace-Beltrami operator A_A4 on M, we can give a 
more intrinsic definition of the spaces H S (A4), by considering 

H 2m (M) = {u£ L 2 (M)/A% t u£ L 2 {M)} , 

which, equipped with the canonical norm 

( 2 -4) \\ U \\ 2 H^(M) = W u \\l 2 (M) + \\^M u \\l 2 (M) ' 

is a Hilbert space. 
We set 

V := {v G H l (M); J M v(x) dM = 0}, 

which is a Hilbert space with the topology endowed by H l {M). 

The condition j M v(x) dM. = is required in order to guarantee the validity of the 
Poincare inequality, 

(2.5) \\f\\l HM) < (Air'UVT/ll!^), for all / G H\M), 

where Ai is the first eigenvalue of the Laplace-Beltrami operator. 

We observe that the problem (|2.2p can be written in the following form 

^ + AU = G(U t ), 
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where 



A 



-/ 
-A M 

is a maximal monotone operator and G(-) represents a locally Lipschitz perturbation. So, 
making use of standard semigroup arguments we have the following result: 

Theorem 2.1. 

• (i) Under the conditions above, problem \2.2\ is well posed in the space VxL 2 (A4), 
i.e. for any initial data {u^u 1 } € V x L 2 (M), there exists a unique weak solution 
of h2.2\) in the class 

(2.6) ueC(R + ;V)nC 1 (R + ;L 2 (M)). 

• (ii) In addition, the velocity term of the solution have the following regularity: 

(2.7) u t G L 2 oc (R + ; L 2 (M)) , 
(consequently, g (m) £ L 2 oc (R+; L 2 (M)) by Assumption \2.1\ 

Furthermore, i/jit ,?/} £ |yni/ 2 (A / () x V} then the solution has the following regu- 
larity 

ueL°° (M+; y n H 2 (M)) n VF 1 ' 00 (M+; V) n VF 2 ' 00 (M+; L 2 (M)) . 



Supposing that u is the unique global weak solution of problem (12. 2h . we define the 
corresponding energy functional by 



dM. 



(2.8) E(t) = \ [ \\u t (x,t)\ 2 + \V T u(x,t)\ 2 

* Jm 1 

For every solution of (|2.2p in the class f|2.6[) the following identity holds 

(2.9) E(t 2 ) - E(ti) = - I ( a(x) g{u t )u t dMdt, for all t 2 > h > 0, 

Jti Jm 

and therefore the energy is a non increasing function of the time variable t. 

3. Main Result 

Before stating our stability result, we will define some needed functions. For this pur- 
pose, we are following the ideas firstly introduced in Lasiecka and Tataru |LA-TA| . For the 
reader's comprehension we will repeat them briefly. Let hbe a concave, strictly increasing 
function, with h (0) = 0, and such that 

(3.1) h (s g(s))) >s 2 + g 2 (s), for \s\ < 1. 

Note that such function can be straightforwardly constructed, given the hypotheses on 
g in Assumption 12.11 With this function, we define 

(3-2) r{.)=h{ — — ). 

meas [hi) 

As r is monotone increasing, then cl+r is invertible for all c > 0. For L a positive constant, 
we set 

(3.3) p(x) = (cl + r)- 1 (Lx) , 
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where the function p is easily seen to be positive, continuous and strictly increasing with 
p(0) = 0. Finally, let 

(3.4) q(x) = x - {I + p)' 1 {x) . 
We can now proceed to state our stability result. 

Theorem 3.1. Assume that Assumptions \2. il and Assumption \2.2\ are in place. Let u be 
the weak solution of the problem \2.2\) . With the energy E(t) defined as in \2. 8\) . there 
then exists a Tq > such that 

(3.5) E(t) < S (Jr ~ l) , Vi>T , 

with lim S(t) = 0, where the contraction semigroup S(t) is the solution of the differential 

t— >oo 

equation 

(3.6) j t S(t)+q(S(t))=0, S(0)=E(0), 

(where q is as given in {3.$ )- Here, the constant L (from definition \3.3\l ) will depend on 

meas(E)(l+||a|| 00 ) ' 



meas(T,), and the constant c(from definition 13. 3\) ) is taken here to be c — ' +/v 



Remark 3.1. If the feedback is linear, e. g., g(s) = s, then, under the same assumptions 
as in Theorem \3.1l we have that the energy of problem \2. 2\) decays exponentially with 
respect to the initial energy, there exist two positive constants C > and k > such that 

(3.7) E(t) < Ce~ kt E(0), t > 0. 

As another example, we can consider g(s) = s p , p > 1 at the origin. Since the function 
p+i 

s 2 is convex for p > 1, then solving 

(3.8) St + S^ = 0, 
we obtain the following polynomial decay rate: 

E(t) < C(E(0)) [£(0)^ +t{p- 1)]=FFT. 

We can find more interesting explicit decay rates in Cavalcanti, Domingos Cavalcanti 
and Lasiecka [CA-DO-LA] . 

4. Proof of Main result 

4.1. Preliminaries. 

We collect, below, some few formulas to be invoked in the sequel. 

Let v be the exterior normal vector field on Ai. For all x € A4, we denote by tt(x) the 
orthogonal projection on the tangent plane T X A4. Any regular vector field q : M 3 — > M 3 
will be split up as follows: 

(4.1) q{x) = q T + (q(x) ■ u(x))u(x), 

where qx = Tr(x)q(x) is the tangential component of q. 
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If tp : R 3 — ► R is a regular function, we have 

(4.2) V<p = d v tpv + V T <f on M, 

(4.3) |V</?| 2 = \d u <p\ 2 + \V T p\ 2 on M, 

where d u , is the normal derivative towards the exterior of A4 and VtV> is t ne tangential 
gradient of 99. 

The Laplace- Beltrami operator Aj^ of a function tp : A4 — ► R of class C 2 is defined by 

(4.4) A_A/fy? := oIwt'Vt^P, 

where divT^Tf, is the divergent of the vector field VtV 9 - 

Assuming that 99 : A4 — > R is a function of class C 1 and q : R 3 — > R 3 be a vector field 
of class C , we have, 

(4.5) / qr ■ VrpdM = — I divqx'pdM, 

Jm Jm 

(4.6) 2<p(qT ■ V^v?) = qr ■ ^ri'-p 2 )- 
From (4.5) and (4.6), we conclude the following formula 

(4.7) 2/ tp(q T -V T tp)dM = / q T • V T (^ 2 ) dM = - / div T qr\<p\ 2 dM. 

Jm Jm Jm 

We observe that in the particular case when m(x) = x — x°, x G R 3 and x° £ R 3 is a 
fixed point in R 3 , we have 

(4.8) V ■ m = 3, diuy my = 2 + (m ■ v)Tr B. 

where B is the second fundamental form of j\A (the shape operator) and Tr is the trace. 
Let tp and m defined as above. We also have, 

(4.9) Vt<P • Vt^t-t • Vt</2 = |Vt<^| 2 + {m ■ z^)(Vt¥? • B ■ Vt</?)- 

The proof of the above formulas can be found in [NE], [LEMlJ, HEM3] and references 
therein. 

Remark 4.1. The sign of B can change in the literature. In our case, we remember that 
B = —dN, where N is the Gauss map with respect to v. 
The formula ( f^.<§| ) can be rewritten by 

(4.10) V-m = 3, div T m T = 2 + 2H(m-v). 
where H = is the mean curvature of j\A . 

We define a continuous linear operator — : H l {A4) — > (i7 1 (A^l)) / , where A4 is a 
nomempty open subset of Ad (sometimes the whole M) such that 



(4.11) 



(-A M tp,^)= [ v T tp-v T ^dM, Vtp^eH^M) 
Jm 
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and, in particular, 



(4.12) (-A M <p, (p )=[ \V T <p\ 2 dM, V<peH\M). 

JM 

The operator — A^+I defines an isomorphism from H l (M) over [H 1 (M)]'. We observe 
that when M is a manifold without boundary, and this is the case, for instance, if M = M, 
we have ^(M) = Hjj(M) and, consequently, [H l (M)\ = H- 1 (M). 

Remark 4.2. It is convenient to observe that all the classical formulas above stated can 
be extended for Sobolev spaces by using of density arguments. 

The proof of Theorem 3.1 proceeds through several steps. 



4.2. An identity. We begin by proving the following proposition 

Proposition 4.2.1. Let M C M 3 be oriented regular compact surface without boundary 
and q a vector field with q = qr + (q • v)v. Then, for every regular solution u of f) 1 . 2 j) we 
have the following identity 

(4.13) / u t q T -V T udM +-[ [ (div T q T ) { |n t | 2 - | V T u\ 2 \ dMdt 
Jm Jo 2 jo J m I ) 

+ [ [ V T u-V T q T -V T udMdt+ f f a(x) g(u t )(q T ■ V T u)dMdt = 0. 
Jo Jm Jo Jm 



Proof. Multiplying the equation of (|1.2p by the multiplier q? ■ Vtu and integrating on 
Mx]0,T[, we obtain 

(4.14) 0= I [ (u tt - A M u + a(x)g(u t ))(q T -V T u)dMdt. 

Jo Jm 

Next, we will estimate some terms on the RHS of identity (|4.14p . Taking (|4.1ip . (|4.6p 
and (|4.7p into account, we obtain 

(4.15) [ [ (-A M u) (q T ■ V T u) dMdt = [ [ V T u ■ V T (qr ■ V r n) dMdt 
Jo Jm Jo Jm 

= [ [ V T u-X7 T qT -V T udMdt + - [ [ q T ■ V T [\V T u\ 2 ]dMdt 
Jo Jm 2 J J m 



l 



JM 2 ,; j M 



T 



Vtu ■ Vrqr • VxudMdt — - I I \Vtu\ 2 divxqr dMdt, 



WAVE ON COMPACT MANIFOLDS 

and, integrating by parts and considering (|4.7I) . we obtain 

i-T 



(4.16) 



(u tt + a(x) g{u t )) (qr • Vx«) dMdt 



o JM 



I ut(qr ■ Vtu) — / ut(qT • Vrut)dMdt 
Jm Jo Jo Jm 



- I / a(x) g (u t ) (qr • Vru)dMdt 
'o Jm 

lT l * 



/ u t {qT-V T u) + - / / (div T qT)\u t \ 2 dMdt 
JM Jo 2 Jo Jm 



4- / / a(x) g (u t ) (q T ■ V T u)dMdt. 

'o Jm 



Combining (I4.14|) , (|4.15|) and ( I4.16P , we deduce ()4.13p , which concludes the proof of 
Proposition 4.2.1. □ 



Employing (|4.13p with q(x) = m(x) = x — x° for some x° E M 3 fixed and taking ((4.8 
and (14. 9j) into account, we infer 



(4.17) 



/ u t m T -V T udM + [ [ \\u t \ 2 - \V T u\ 2 \ dMdt 
Jm J o Jo Jm l j 

+ / f [|V T u| 2 + (m- i/)(V T «- B -V T u)]dMdt 
Jo Jm 



-}- / / a(x) g(u t )(m T ■ V T u)dMdt = 0. 
'o Jm 



We have the following identity: 

Lemma 4.2.3. Let u be a weak solution to problem (|1.2p and £ € C l {M). Then 



(4.18) 



ut^u dM 



M 



[ £\u t \ 2 dMdt- [ [ £\V T u\ 2 dMdt 
o Jm Jo Jm 

[ (VTU-V T £)udMdt- [ [ a(x) g(u t )£udMdt. 
o Jm Jo Jm 



Proof: Multiplying the first equation of (jl.2p by £it and integrating by parts we obtain 
the desired. □ 
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Substituting £ = \ in (j4.18j) and combining the obtained result with identity (14.170 we 
deduce 



(4.19) 



/ 

JM 



Ut rriT • Vtu dM 



+ 



1 



n 



ut u dM 



M 



o 



E(t)dt+ [ [ a(x) g(u t )(m T -V T u)dMdt 
JO JM 



1 

+ 2 



T 



a(x) g{ut)udMdt 



o JM 

T 



-J J ( m - v) H {\ u A 2 - \V T u\ 2 }dMdt. 

[ [ (m- v)(y T u- B -V T u) dMdt. 
JO J M 



Analysis of the terms which involve the shape operator B 

Let us focus our attention on the shape operator B : T X M — > T X M. There exist an 
orthonormal basis {ei,e2} of T X M such that Be\ = k\e\ and Be2 = &2 e 2- k% an d &2 are 
the principal curvatures of M at x. The matrix of B with respect to the basis {ei,e2} is 
given by 



B 



ki 
k 2 



Setting Vt« = (£,??) the coordinates of Vt^ in the basis {e\,e2}, for each x € M, we 
deduce that 



(4.20) 

Then, from (|4.20p . we infer 



Vtu ■ B ■ Vtu = k\£ 2 + k2rf. 



(4.21) 



(m • 
= (m • z/) 



(V T u • S • V T u) - -Tr(B)\V T u\ 2 



(h - k 2 ) 2 . (&2 - fel) 2 
o ? o ^ 



Remark 4.3. Observe that this is the precise moment that the intrinsic properties of the 
manifold M appear, that is, we strongly need that the term — Jq Jj^( m ■ v)Hul dM dt 
lies in a region where the damping term is effective. Remember that the damping term is 
effective on an open set M* which contains M\ U* =1 Moi- So, assuming that H < and 
since m(x) ■ v{x) < on Mq, we have 



[ [ (m ■ u)H \u t \ 2 dMdt < 0. 
Jo J Mo 
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In addition, supposing that Moi is umbilical for every i = 1, . . . , k, then, having \4-21 ) 
in mind, we also have that 

[ [ {m-u)[H\V T u\ 2 -(V T u-B -V T u)] dMdt = 0, i = l,...,k. 

JO JMoi 

More generally, assuming that the principal curvatures k\ and k 2 satisfy \k\(x)—k2{x)\ < 
£i (here, £j is assumed sufficiently small) for all x £ Mqi, i = 1, • • • , k, we deduce that 

k ,. T . 

V / / (m • v) [H\V T u\ 2 - (V T n • B ■ V T u)] dMdt 

■ =1 Jo J Mo, 
k „ T „ 

<y^/ / \(m-v)\\k 1 -k 2 \\f + n 2 \dMdt 
i=1 Jo J Mm 

k „ T . k .. T 

<y2Ri£i / \V T u\ 2 dMdt <2Vi?i£i/ E(t)dt, 

i=1 Jo Jm m i=1 ^0 

w/iere i?j = max xg: ^-||x - a; || K 3. 

Set 7W 2 = M\ Uf =1 A^oi- In the case where .Moi are umbilical, recalling (14.19P taking 
(|4.2ip and Remark 4.3 into consideration, we deduce 

T 1 r r i T 

-^udM 



(4.22) I E(t)dt < - f u t m T -V T i 
Jo Um 



- T 1 






/ ut u dM 


-o 2 


J M 



+ 



1 

+ 2 



/ / {m-v) [H\V T u\ 2 - (V T u ■ B ■ V T u)] dMdt 

Jo Jm 2 

f f (m ■ v)H \u t \ 2 dMdt 

Jo Jm 2 

a(x) g(ut)(m,T ■ Vxu)dMdt 
a(x) g(ut)udMdt. 



T 

JM 



T 

JM 



In the general case, the unique difference in the proof is that the term f Q E(t) dt that 
appears on the LHS of (|4.22p will be multiplied by a positive constant C, provided that 
we consider £j small enough. For simplicity we shall assume that C = 1. 

We will denote 



(4.23) 



X 



ut rriT ■ Vtu dM 



M 







r 


/ ut u dM 






JM 






Next we will estimate some terms in (|4.22p . Let us denote: 

(4.24) R : = max [|m(a;)[|Rn = max ||x — x ||rh. 

x^M x&M 



Estimate for 1% := J Q T J M a{x) g{ut){rriT ■ Vxu)dMdt. 
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By Cauchy-Schwarz inequality, taking (|4.24|> into account and considering the inequality 
oh < |jj + r]b 2 , where r\ is a positive number, we obtain 

(4.25) |/i| < IM^ZlMj / /" a (x)\g(u t )\ 2 dMdt + 2t] [ E{t)dt. 

V Jo Jm Jo 

Estimate for I2 = \ Jq J m a(x) g(ut)u dAidt. 
Similarly we infer 

Hall A -1 f T r r T 

(4.26) \I 2 \ < UL ™[ M) 1 / / a{x)\g{u t )\ 2 dMdt + 2r } / E{t)dt, 

16r? J JAf Jo 



where Ai comes from the Poincare inequality given in (j2.5f) . 

Choosing r/ = 1/8 and inserting P~23j) . P~25j) and (|4~26|) into ([4T22]) yields 

(4.27) -/ E(t)dt < |x| + Ci/ / a(i)( 9 (u ( )) 2 ^rft 

^ Jo Jo Jai 

+ Cl / / [|V T u| 2 + a(x)u 2 ]dMcft 
JO JA4 2 



where 



Ci := maxlllallioo^)^"^^ 1 + 8i? 2 ], + \H\R, RlH^ 1 } , 



||-B|| = sup \B X \, with (S^l = sup l-B^I- 
x£M {v£T x M;\v\=l} 

It remains to estimate the quantity f M ^ \Vtu\ 2 dMdt in terms of the damping term 

lo -I'm [ a ( x ) \9( u t ) I 2 + a ( x ) I - "* I 2 ] dAidt. For this purpose we have to built a "cut-off" func- 
tion r] e on a specific neighborhood of M.2- First of all, define fj : R — ► R such that 

( 1 if x < 

77(0;) = I (x - l) 2 if x G [1/2,1] 
( if x > 1 

and it is defined on (0, 1/2) in such a way that fj is a non-decreasing function of class C . 
For e > 0, set %(x) := fj{x/e). It is straightforward that there exist a constant M which 
does not depend on e such that 

\v'e(x)\ 2 < M 

fj £ (x) ~ e 2 

for every x < e. 

Now let £ > such that 



i=l 



:= {x G AI; dist(x, |J clMoi) < e} 

ghborhood oi 

r/ e : M — > R as 



is a tubular neighborhood of Ui=i ^JWoi and u; £ := Cj £ U AI2 is contained in hi*. Define 



1 if x G Al 2 

r? £ (x) = { fj £ (d(x,M 2 )) if xGw £ \-^2 
otherwise. 
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It is straightforward that i] £ is a function of class C 1 on M due to the smoothness of 
dM 2 and dio £ . Notice also that 



(4.28) 



iVrTfeQzQI 2 _ \fj' £ (d(x,M 2 ))\ 2 < M 
r] £ (x) fj e (d(x,M 2 )) ~ £ 2 



for every x E u e . In particular, £ L°°{uj £ ). 

Taking £ = rj £ in the identity (|4.18p we obtain 
rT 



(4.29) 



Jw E 



rj £ \V T u\ 2 dMdt 



1 T 



utur] £ dM 



+ 

JO Juj, 



T 



Jul 



u{Vtu ■ VtVs) dMdt 



i] £ \u t \ 2 dM 

-fj 

JO Juj, 



a(x) g(ut)urj £ dMdt. 



Next we will estimate terms on the RHS of (|4.29|) . 
Estimate for K \ := Jq ri £ \ut\ 2 dMdt 

From (|2.3|) . since r\ £ < 1 and uj £ C M*, where the damping lies, we deduce 



(4.30) 



a(x) u t dM dt. 



JM 



Estimate for K 2 := — J w a(x) g(ut)urj £ dMdt. 

The Cauchy-Schwarz inequality, the inequality ab < j^a 2 + ab 2 and (|2.5p yield 



(4.31) \K 2 \ < 



A i ll a lk o °(M) 



4a 



T 

JM 



a{x)\g{u t )\ 2 dM + 2a / E(t)dt, 
Jo 



where a is a positive constant. 

Estimate for K3 := J Q T u(Vtu ■ VTf]e)dMdt. 

Considering (14. 28ft and applying Cauchy-Schwarz inequality, we can write 



(4.32) \K 3 \ < \f T \f r l£ \V T u\ 2 dM+ I' 

^ JO Uuj e Jw t 



^l\u\ 2 dM 



dt 



< 



■q £ \V T u\ 2 dM + 



M 



\u\ 2 dM 



dt. 



Combining f|4.29|) - (|4.32|) we arrive to the following inequality 



(4-33) - 



/ r] £ \V T u\ 2 dMdt < \y\ + - 

Jw E 



\L°°(M) 



4a 



T 



+ 



a(x) \g{u t )\ 2 dM 

t x M I' 1 f 

2a E{t)dt + — / \u\ 2 dMdt, 
Jo 2e ,/ J u 



JM 
T 



-1 



T 

JM 



a(x) u t dM dt. 



14 M. M. CAVALCANTI, V. N. DOMINGOS CAVALCANTI, R. FUKUOKA, AND J. A. SORIANO 
where 



(4.34) 



u t urj £ dM 



Thus, combining (j4.33j) and (I4.27|) . have in mind that 



If [ \V T u\ 2 dMdt<\ [ [ V£ \V T 

1 JO J MO 1 JO JLUr 



-u\ dMdt 



io JM 2 

and choosing a = 1/1 6Ci we deduce 



(4.35) ±J E{t)dt<\ X \ + 2C 1 \y\ 

+max{Ci,8Cf\^ 1 \\a\\ L oo^ M ),2Cia 1 } [ [ [a(x) \g(u t )\ 2 + a(x) \u t \ 2 } dMdt 

Jo J M 



MCi 



r 



Ju) 



Id dM dt. 



On the other hand, from (|4.23p . (|4.34|) and (|2.9p the following estimate holds 
(4.36) |xl+2C 2 |^| < C(E(p)+E(T)) 

= C 2E(T)+ [ I a{x)g{u t )u t dM , 
jo Jm 

where C is a positive constant which depends also on R. 
Then, (|435j) and (j436j) yield 



(4.37) TE(T) < [ E(t) 

Jo 



dt 



< CE(T) + C 



[a(x) \g(u t )\ 2 + a(x) \u t \ 2 } dMdt 



JM 



+ C 



Jlu 



\u\ dM dt, 



where C is a positive constant which depends on ao, Ai,i?, \H\, \ \B\\ and 

Our aim is to estimate the last term on the RHS of (|4.37p . In order to do this let us 
consider the following lemma, where To is a positive constant which is sufficiently large 
for our purpose. 



Lemma 4.1. Under the hypothesis of Theorem \3.1\ and for all T > Tq, there exists a 
positive constant C(Tq, E(0)) such that if(u,u t ) is the solution of jlLOfy with weak initial 
data, we have 

(A.38)£ \\u(t)\\ 2 LHM) dt < C(T ,E(0)) j^(a(x)g 2 (u t )) + a(x)u 2 ) dMdt} . 
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Proof: We argue by contradiction. For simplicity we shall denote u' := ut- Let us 
suppose that (14. 38ft is not verified and let (-Ufc(O), u' k (0)} be a sequence of initial data where 
the corresponding solutions {/Ufc}fceN of (11. 2} with E k (0), assumed uniformly bounded in 
k, verifies 



So \\vk(t)\\h(M) 



dt 



(4.39) lim -= = +oo, 

^ + °° Jo S M H x ) 9 2 K) + <x) u'Z) alM dt 

that is 

(440) lim So Sm W 9 2 K) + a{x) u' 2 ) dMdt _ q 

k ^+°° So WMt)\\h {M ) dt 

Since E k (t) < E k (0) < L, where L is a positive constant, we obtain a subsequence, still 
denoted by {u k } from now on, which verifies the convergence: 

(4.41) u k -± u weakly in il^Er), 

(4.42) Mfc — u weak star in L°°(0, T; V), 

(4.43) < u' weak star in L°°(0,T; L 2 (M)). 

Employing compactness results we also deduce that 

(4.44) u k -> it strongly in L 2 (0, T; L 2 (7W)). 

At this point we will divide our proof into two cases, namely, a / and u = 0. 

(i) Case (I): u ^ 0. 

We also observe that from (|4.40|) and (|4.44p we have 

(4.45) lim / f (a(x) g 2 (u' k ) + a{x) u' k 2 ) dMdt = 

fc^+oo J Q J M 

Passing to the limit in the equation, when k — > +oo, we get, 

' u# - A M u= on A4 x (0, T) 
u t = 0on^ t x(0,T), 

and for ut = v, we obtain, in the distributional sense 

(v tt -A M v= Oon M x (0,T), 
| u= on M* x (0,T). 

From standard uniqueness results we conclude that v = 0, that is, Ut = Returning to 
(|4.46p we obtain the following elliptic equation for a.e. t £ (0, T) given by 

A_A/f u = on A4 
ut = on M, 

which implies that u = 0, which is a contradiction. 

(ii) Case (II): u = 0. 



(4.46) 
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Defining 



Cfc := 



\u k \ 2 dM dt 



JM 



1/2 



(4.47) 

and 
(4.48) 
we obtain 

(4.49) J T f M \u k \ 2 dM dt = f M ^fdM dt = jz ft f M \u k \ 2 dM dt = 1. 

k k 

Setting 

E k (t) 



u k := — u k , 
c k 



1 



M 



\vf k \' z dM + | / [V?Z fc | 2 d.M, 



we deduce automatically that 
(4.50) 



E k (t) 



E k (t) 



Recalling (|4.37p we obtain, for T large enough that 



E(T) < C 



JM 



(a(x) g 2 (u t ) + a{x) u 2 t ) dM dt + 



\u\ dMdt 



o JM 



and employing the identity E{T) — E(0) = — J M a(x) g(ut) u t dM. dt, we can write 



E(t) < E(0) < C 



(a(x) g 2 {u t ) + a{x) uf) dM dt + 



o JM 



I til dMdt 



.7.M 



for all t G (0, T), with T large enough. The last inequality and (|4.50p give us 



(4.51) E k (t) :=^<C 



£f M (a(x)g 2 (u> k ) + a(x)u> 2 ) 
lo I M \u k \ 2 dMdt 



+ 1 



From (|4.40p and (|4.51h we conclude that there exists a positive constant M such that 



E k (t) :-- 



E k (t) 



< M, for all t G [0, T] and for all k G N, 



that is, 



(4.524 / \vf k \ 2 dM + - \Vu k \* dM < M, for all t G [0, T] and for all k G N. 

For a subsequence {ufc}, we obtain 

(4.53) u k weak star in L°°(0, T; V), 

(4.54) u' k u' weak star in L°°(0, T; L 2 (.M)), 

(4.55) u fc — > u strongly in L 2 (0,T; L 2 (M)). 
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(4.57) 



We observe that from (|4.40p we deduce 
(4.56)im [ T f a ( g ) g2 K) dMdt = and lim ft a(x) \u' k \ 2 dM dt = 0. 

k^+ooJ Q J M C k k^+ooJ Q J M 

In addition u k satisfies the equation 

vfl - A M u k + a{x) =0 on M x (0, T). 

Ck 

Passing to the limit when k — ► +oo taking the above convergence into account, we 
obtain 

( u"-k M u= onA4x(0,T), 
vf = onM*x(0,T). 

Then, v = Ut verifies, in the distributional sense 

vu — Ayvi v = on M. 

v = on M*. 

Applying uniqueness standard results it results that v = Tit = 0. Returning to (|4.57p 
we obtain, for a.e. t G (0, T) that 

Am u = on M 
ut = on M., 

from what we deduce that u = 0, which is a contradiction in view of (|4.49p and (|4.55p . 
The lemma is proved. □ 

Inequalities (|4.37p and (|4.38p lead us to the following result. 

Proposition 5.2.2: For T > large enough, the solution [u, ut] of (|2.2p satisfies 



r 



(4.58) E(T)<C / / a{x)\u t \' + a{x)\g{u t )Y 

Jo J Mi 



dMdt 



where the constant C = C(Tq, E(0), oq, Ai, R, \\B\ 



7?) 



4.3. Conclusion of Theorem 3.1. In what follows we will proceed exactly as in Lasiecka 
and Tataru's work |LA-TA] (see Lemma 3.2 and Lemma 3.3 of the referred paper) adapted 
o our context. Let 

S a = {(i, x) G Si/ \ut\ > 1 a. e.}, 
S/3 = Ei\E a . 
Then using hypothesis (Hi) in Assumption 12.11 we obtain 

(4.59) f a(x) (\g (u t )] 2 + (n t ) 2 ) dS Q < (A; -1 + K ) J a(x)g (u t ) u t dZ a . 

Moreover, from (|3.ip 



(4.60) / a{x)([g{u t )] 2 + {u t ) 2 )dY 1 p<{l + \\a\\ 00 ) h(a(x)g (u t ) u t ) d^p 
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Then by Jensen's inequality 

(1 + Halloo) / h(g(ut)ut)dEp < (1 + ||a||oo)meas (£) /i ( 3_— / a(x)g (u t ) u t dE J 

Jn \meas(Z)Jx ) 

(4.61) = (1 + \\a\\oo)meas (T,) r (J a(x)g (u t ) u t dJ?j , 

where r (s) = h ( meffl g( S ) ) is defined in (|3.2[) . Thus 

y o(ar) ([g(u t )) 2 + (u t f)dZ < (k~ l + K) J a(x)g (u t ) t cEi 
(4.62) +(1 + 110(100)771605 (S) r a(x)g (ut) UfdT^j . 



Splicing, together, (|4.58[) and (|4.62p . we have 

'l + llol 



00J JT, 



E{T) < (l + ||a||oo)C 

(4.63) +meas(T,)r^J a(x) g (ut) utdY^J 

where Kq = k~ l + K. Setting 

C meas (E) (1 + ||a||oo) ' 
M 



g (u t ) utdZi 



meas (£) (1 + ||a||oo) 

we obtain 

(4.64) P[E(T)\ < I a(x)g( Ut )u t dZ = E(Q)-E(T), 

Jt, 



where the function p is as defined in (|3.3f) . To finish the proof of Theorem 3.1, we invoke 
the following result from I. Lasiecka et al. |L A-TAj : 

Lemma B: Let p be a positive, increasing function such that p(0) = 0. Since p is 
increasing we can define an increasing function q, q(x) = x — (I + p)" 1 (x) . Consider a 
sequence s n of positive numbers which satisfies 

+ p{s 

m+l) — ^in- 
Then s rn < S(m), where S(t) is a solution of the differential equation 

±S(t)+q(S(t))=0, S(0)=s . 
dt 

Moreover, if p{x) > for x > 0, then lim S(t) = 0. 

With this result in mind, we replace T (resp. 0) in (|4.64p with m(T + 1) (resp. mT) 
to obtain 

(4.65) E(m(T + 1)) + p (E(m(T + 1))) < E(mT), for m = 0, 1, .... 
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Applying Lemma B with s 



E(mT) thus results in 



(4.66) 



E{mT) < S(m) 



m = 0, 1, .... 



Finally, using the dissipativity of E(t) inherent in the relation (|2.9I) , we have for t = 
mT + t, < r <T, 



where we have used above the fact that S(.) is dissipative. The proof of Theorem 3.1 is 
now completed. 

4.4. Further Remarks. From the above procedure, we can construct a wide assortment 
of compact surfaces by jointing pieces of different kind of surfaces. However, according 
to the position of the observer (point xq) the dissipative and the non dissipative areas 
can change drastically. To illustrate this, let us consider the Catenoid or the Trinoid (see 
figures 2 and 3 below) that are minimal surfaces, that is H = 0. 



Considering a strategic piece of one of these surfaces we can construct another compact 
surface according to the figure 5 above. Remember that the non dissipative regions must 
occur where m(x) ■ v{x) < 0, H < and simultaneously we are forced to consider \k\ — 
sufficiently small (by parts) . The dissipative area must contain strictly the closure of the 
points x € M such that m(x) ■ v(x) > 0. It is not difficult to see that the non dissipative 
area in the figure A can be located near the top and/or near xo while the non dissipative 
area in the figure B can be located in the middle in the middle of surface and/or near xo, 
assuming evidently that k\ ki on these non dissipative areas. 



(4.67) 




for t > T 




Figure 2: Catenoid 
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Figure 4: Trinoid 




x 

Figure 3. 
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